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APPLIED SCIENCES 


GLOTHERMAL FIELD STRATIGRAPHIC PRESSURE GRADIENTS ANALYZED 


Beijing TANKUANG GONGCHENG [PROSPECTING ENGINEERING) in Chinese No 2, 22 Apr 80 
pp 1-8 


[Article by Wang Zhongying [3076 0112 5391), Zhou Shaoguang [0719 7300 0342) and 
Tang Songran [3282 2646 3544): "Categorization of Stratigraphic Pressure Gradi- 
ents in Geothermal Fields and the Laws of the Stratigraphic Gradient in the Yang- 
bajing Geothermal Field") 


[Text] Abstract: Gecthermal reservoirs, whose hot fluid can be 
continuously replenished, are open-type structures, in contrast 
to ofl and natural gas reservoirs which are basically of the 
closed type. Accordingly, the stratigraphic pressure classifi- 
cation methods used in drilling ofl and gas wells are basically 
not applicable to geothermal boreholes. This article attempts 
to classify stratigraphic pressures in geothermal wells and 
studies the geothermal gradient laws of the Yangbajing Geother- 
mal Pield. The equations presented have been tested in 3 years’ 
drilling experience at the Yangbajing field against the actual 
date from more than 10 wells. 


1. Statement of the Problem 


The Yangbajing Geothermal Field, Tibet, is our country's first high-temperature 
moist steam geothermal field. Because the reservoir is shallowly buried (already 
encountered at 20-30 meters) and is in a thick stratum (several hundred to more 
than 1,009 meters), while the cover layer is thin and brittle, it is very diffi- 
cult to control the pressure within the well. If the specific gravity of the mud 
is too low, it may result in gushing, while if it is too high it may result in 
infiltration followed by gushing. An even more difficult problem is that if the 
uncased bottom section of the well is a bit too long, it may result in gushing in 
the upper section and infiltration below. In the first 2 years of prospecting, 
serious gushing caused stratigraphic subsidence, well-head explosions or surface 
cracking, resulting in the abandonment of 4 wells during the drilling process, 
drilling work was unable to continue, since there was 4 grave threat to construc- 
tion safety. 


The experience of well drilling at the Yangbajing Geothermal Field has posed a 
series of problems: 














A. Wow great 4 pressure is involved in serious blowoute’ Can they be controlled 
with mud? 


5B. On the basis of foreign experience with geothermal wells and well suppression, 
a method of pressurized cold water control of blowoute has been published; why 
can thie be successfully used on some wells at the Yangbajing field but not on 
others’ How should the well pressurizing method be chosen at the Yangbajing 
field? 


C. Why cannot the length of the uncased bottom section be too long, particu- 
larly in a shallow section, without the danger of gushing above and infiltration 
below. How should the position of the casing point be chosen’ 


D. How should the mud specific gravity be chosen in order to control the pres- 
sure inside the wel! and prevent both blowouts and infiltration? 


E. Why do some wells gush naturally after completion, while others must be sub- 
jected to gas lift or some other special method to induce gushing? Why is the 
pressurized cold water method of inducing flow which is used extensively abroad 
not applicable to some wells at the Yangbajing field? 


Although these engineering probleme differ from each other, they are all con- 
strained by a common factor, namely stratigraphic pressure and the pressure 
gradient. 


The problem of stratigraphic pressure and pressure gradients is a serious one 
which urgently requires solution. It has long received serious attention in the 
petroleum and natural gas well drilling world, ae many methods for prediction 
before or during drilling have been developed,’ such as the seismic method, which 
uses data on the propagation of manmade seismic waves, the sonic logging and 
inductance logging methods based on logging data from wells with known character- 
istics, the d-coefficient method which uses engineering data from the drilling 

of the well, and the like. In addition, on the basis of the conditions of 
formation and accumulation of ofl and gas reservoirs, strata have been classi- 
fied as high-pressure, normal-pressure and low-pressure types. The methods and 
the classification have reference value for the drilling of oil and gas wells, 
but because the stratification, heat reservoir characteristics and reservoir 
structure of geothermal fields differ widely fron those of oil and gas strata 
and reservoirs, oil and gas well stratigraphic pressure classification and pre- 
diction methods are basically inapplicable to geothermal wells. To date no one 
in this country or abroad has attempted a classification or caloulegion of strati- 
graphic oressure in geothermal wells. We have published an article’ providing a 
simple introduction to stratigraphic pressure gradients in geothermal fields as 
one of the problems of geothermal well drilling technology. The present article 
will supplement and extend the previous one. Because of limited space, certain 
mathematical manipulations and universally known concepts are omitted. Inter- 
ested readers may consult another article in the collection from the Second 
Technical Conference of the Prospecting Engineering Committee of the China Geo- 
logical Society.’ 

















2. Classification of Stratigraphic Pressure Gradients 


In order to clarify the pressure-bearing states of strata in a way that is help- 
ful in choosing the proper engineering approach, we must analyze the data in 

terms of pressure. Because stratigraphic pressure gradients are closely related 
to engineering questions, *o the sedimentary environments of strata and to later 
compaction or pressure release effects, the stratigraphic pressure classifications 
are all based on the magnitude of the stratigraphic pressure gradient. 


In of] or gas wells, strata are generally classified as high, normal or low 
pressure types on the »asis of whether or not there have been subsequent com- 
pacticn or pressure leakage effects. 


In high pressure strata, the pressure gradient is from more than 0.107 to 0.23 
kg/cm‘ =m; 


in normal-pressure strata the pressure gradient is from 0.10 to 0.107 kg/cm*-m; 
in low-pressure strata the pressure gradient is from 0.09? to 0.055 kg/cm? =m. 


With the original sedimentary state unchanged, strata which have not been tightly 
sealed and compacted or had pressure leakage are normal-pressure strata; because 
the pressure gradient of fresh water is 0.10 kg/cm*-m and that of sea water is 
0.107 kg/cm*=-m, in fresh water facies (lacustrine) or sea water facies (marine) 
sedimentary strata the pressure gradient of a normal-pressure stratum should be 
between the gradients of fresh water and sea water, i.e., between 0.10 and 0.107 
kg/cm? =m. 


Stratigraphic pressure gradients higher or lower than those of normal-pressure 
Strata are the result of subsequent sealing and compaction or pressure leakage 
and migration. The upper limit of the stratigraphic gradient in high-pressure 
strata, which is 0.23, is determined from the cap rock pressure gradient. For 
loose, unconsolidated strata, the cap rock pressure gradient is equivalent to 

the maximum internal pressure that can be borne by the stratum. If it exceeds 
0.23, there is a danger of forcing through the cap rock. This is why 0.23 is 
defined as the upper limit; the numerical value is obtained by calculation. 





The factors and strata producing thermal reservoirs are different from those 
producing oil and gas reservoirs. Geothermal reservoirs are usually associated 
with volcanic rock, and the reservoir structure is generally of an open type. 
The cause of formation of hot fluid is that cold surface water penetrates under- 
ground through cracks, and after heat exchange with deeply buried hot rock bodies, 
it rises again through structural cracks to some shallow ly-buried porous or 
fissured stratum, where part may still pass to the surface and be lost. This 
type of system is an open system, and there is no question of later compaction. 
Accordingly, the stratigraphic pressure classification method based on subse- 
quent compaction effects which is used for oil wells has no theoretical or real 
meaning for geothermal wells. 








However, geothermal reservoir strata actually do have different pressure states. 
The deciding fac. or is not subsequent compaction, but rather the temperature of 
the reservoir and the depth at which it is located, High pressures can occur 
particularly in shallowly-buried reservoirs. More deeply buried ones are 
generally of normal or low pressure. 


On the basis of blowout control and infiltration prevention characteristics of 
geothermal wells, we can make the classification of geothermal well strata 
shown in Table 1. 

Table 1. Stratigraphic Pressure Classification of Geothermal Strata 


Stratigraphic Stratigraphic 











pressure pressure gradient 

category (kg/cm* =m) Well construction characteristics 

High pressure Over 0.12 to 0.23 Heavy mud; blowout prevention equipment; 
easy spontaneous gushing 

Normal 0.10 to 0.12 Plain water or nonheavy mud; blowout 

pressure prevention equipment; easy spontaneous 
gushing 

Low pressure More than 0, Foam or air, or use water with alternate 

less than 0. infiltration and drilling 

No pressure Zero or less Foam or air, or water with alternate 
infiltration and drilling; permanent 
stoppage 





——_— —— « —— ———— ee ee 





When the pressure balancing well-drilling method is used, Table 2-1 can be used 
for suitable selection of drilling fluids, including the following series of 
drilling fluids of different specific gravities: heavy mud + light mud + plain 
water + foam (air dispersed in water) + fog (water vapor dispersed in air) + 
air. These fluids can be used in the range between high-pressure strata with 
gradients of 0.23 down to zero-pressure strata with gradients of zero. If we 
use powdered galena as a specific gravity increaser, the fluid can even be used 
in super-high-pressure strata with pressure gradients as high as 0.33-0.36. But 
such strata have not yet been found in geothermal fields, and are quite rare even 
in ofl and gas wells. 


3. The t(Z) Function Relating Stratum Temperature and Depth 


We discuss below the fact that in hot water depth-averaged specific gravity 
calculations and stratigraphic pressure and pressure gradient calculations we 
mus’ understand the law governing the variation of stratum temperature with 
depth; we call this law the function t(Z): 


t = t(Z) (3-1) 


where Z is the depth and t is the stratum temperature. The derivative dt/dZ is 
the stratum temperature gradient. 











Hecaune of the complexity of subterranean conditions, the function t(Z) has many 
forms. Here, on the basis of moist-steam type geothermal fields we can make an 
analysis of possible forms of the function ¢(Z) in terme of heat transfer theory 
and propose some approximate calculation methods. 


According to heat transfer theory, heat transfer is limited to three types: 
radiation, convection and conduction. Each type has its own particular condi- 
tions. Radiative heat transfer requires a relatively high heat source tempera- 
ture to reach any substantial level, while thermal reservoirs are generally at 
rather low temperatures, usually not exceeding 200-300°, so that radiative trans- 
fer can be ignored in comparison with the other two types of heat transfer. 
Convective heat transfer requires the presence of a fluid, and the transfer pro- 
cess is accompanied by mass displacement. Accordingly, this type of heat trans- 
fer exists only in water-containing strata or heat reservoir strata. Conductive 
heat transfer relies on the thermal motion of molecules of solid matter, and it 
does not require transfer of mass. Accordingly, theoretically this type of heat 
transfer should exiet in any type of stratum, but because rocks and strata are 
poor conductors of heat, the heat transfer coefficient is very low, so that in 
water-cont.ining strata the effect of conduction can be neglected in comparison 
with that of convection, while in nonwater-containing strata or cover layers, 
because the effect of convection is very small, the effect of conduction cannot 
be neglected, and it is the main heat transfer type in nonwater-containing strata. 





Thus, on the basis of the heat transfer characteristics of strata we can very 
roughly classify geothermal field strata as convective or conductive. The former 
type are water-containing layers, such as heat reservoirs and cold-water~containing 
strata. The latter are nonwater-containing strata such as the cap rock, water- 
enclosing strata of thermal reservoirs and other nonwater-containing strata. 

Below we discuss the thermal gradient characteristics and t(Z) function for 

these two types of strata. 


A. Conductive Strata 
The decline of geothermal flow with time is very slow, and accordingly for heat 


conduction we can use the Fourier equation for stable heat flow: 


\ 
q"* 7-7, Zo (t - to) (3-2) 


Rewriting this formula and letting 


» S - 
K : (3-3) 


we obtain the formula for t(Z): 
t = to +K (Z = 29) (3-4) 


where q is the quantity of heat (kca!/m“-hr) flowing through 1 square meter of 
surface in an hour. Because this is stable heat flow, and relative to the well 


wi 








segmenc 2 = ‘9 for which the calculation is being performed the heat reservoir 
area is very large, so that the area through which heat flows between 7,9 and 7 
can be considered approximately constant, we can assume that q is approximately 
constant; A ie the conductivity of the stratum (kcal/mhr:deg), and for strata 
at the same level we can consider it to be constant. K is the stratigraphic 
thermal gradient (°C/meter), tg and t are the temperatures (°C) of a known 
point and the calculation point; Zo and Z are the depth of the known point and 
the calculation point (meters). 


For a specific stratum we can consider q and \ to be approximately constant, so 
that K is also a constant. Accordingly formula (3-1l)a is a straight line 
function, 


Because the conductivity A of rock is very small and the unit heat flux q is 
much larger than for an orlinary stratun, the temperature gradient K of the cap 
rock is rather large. A high temperature gradient is a characteristic of cap 
rock. 


B. Convective Strata 


The pores and fissures in convective straia are full of fluid. In general, be- 
cause of the natural convection of hot water, the upper and lower parts of a 
given stratum have rather uniform temperature, which is to say that the tempera- 
ture gradient K = dt/dZ approaches zero. This is the general characteristic of 
convective strata. 


The function t(2) is a straight line function: 
t = const. (3-1)b 


However, when there is entering cold water, the temperature law for the transi- 
tion zone in which cold water and hot water are in contact is relatively com- 
plex. It is affected not only by the position of entry of the cold water and 
the runoff flow rate, but also by the degree of fissuring of the stratum, and 
it is difficult to propose a precise and general formula: we can only use 
approximate methods. We used a trapezoidal approximation approach, using the 
straight-line function (3-l)a to give approximately the function t(Z) for the 
transitional zone. This may involve a rather large error in the transition 
zone, but because its thickness is much less than the depth of the well, the 
error can be neglected in engineering calculations regarding depth-averaged 
specific gravity and stratigraphic temperature gradient. 


To sum up the foregoing, for both conductive and convective layers the func- 
tion t(Z) of a s mple stratum can be expressed by the formula 


t = to + K (Z = Zo) (3-l)a 


For conductive cap rock and convective transitional zones, K will be rather 
large and will be positive for temperature increasing with depth and negative 
for decreasing temperature; for most reservoirs K is close to zero. 


























Figure 1. The different characteristics of the function t(Z) and the pressure 
gradient K indicate that: reservoir C; is affected by cold water runoff and 
is a transition zone; collector Cy is not affected by cold water runoff or is 
only slightly affected; C, and Cy are both gravel reservoirs; reservoir B is 
a silica sandstone reservoir; and reservoir A is a deep fissured reservoir. 





Key: 

a. ca’ egory k. hot water 

b. typical stratum 1. conductive 

c. zone number m. water-impermeable stratum 
d, temperature, t(Z) n. convective 

e. conductive o. Reservoir Co 

f. cover layer p. relatively large 
g. large q. smail 

h. convective (transition zone) r. convective 

i. Reservoir C, s. Reservoir B 

j. cold water t. Reservoir A 


Figure 1 shows the characteristics of the function t(Z) and the thermal gradi- 
ent for typical etrata. The graph of these characteristics is useful when using 
well temperature data to determine the true temperature of the heat reservoir. 
Because the temperature curve for the well interior reflects the temperature of 
the fluid in the well (and not that in the stratum), if no analysis is made and 
the well temperature is taken as the reservoir temperature, the error may be 
very large. If in the reservoir section of the well effects are produced by 
cyclic cooling by the mud and vertical heat convection, the temperature inside 
the well or the temperature of the stratum near the well wall will be lower 

than the true temperature of the heat reservoir; moreover, in a well section in 
a layer impermeable to water or in a cover stratum, owing to gushing or vertical 
heat convection, the temperature inside the well may be much greater than the 
true temperature of the stratum itself, a phenomenon which is especially apparent 





in th. section near the reservoir. Accordingly, in order to minimize the differ- 
ence between temperature measurement data and the reservoir t «perature, before 
taking temperature readings there must be a thorough purging (generally at ieast 
24 hours), and the temperature must be read immediately after closing off the 
well. Moreover, for fissured reservoirs, temperature data in the fissured srea 
is the best indication of the true reservoir temperature. 


4. The Relationship Between the Specific Gravity and Temperature o! Water 


The formula expressing the telationsh >» octween the specific gravity and the 
temperature of pure water can be obt. sed from saturated steam tables.” For 
convenience in calculating stractigra, ovperature gradient formulas, we use 
an approximate calculating formola: 


r' = 1.020 - 0.040° 1072t = 0.020 - 107%? (4-1) 


where t is the temperature (°C) and ty is the specific gravity of pure water. 


The relative error « of approximate formula (4-1) is: 


not greater than 1.10 percent when t = 40-300°C; 
not greater than 0.6 percent when t = 10-280°C; 
not greater than 0.35 percent when t = 80-200°C. 


[hese errors have no effect at all on engineering calculations of the strati- 
graphic pressure gradient. 


A real geothermal fluid will contain large quantities of soluble minerals and 
suspended microscopic particles, and accordingly the specific gravity will be 
greater than for pure water. In the thermal fluid of the Yangbajing Geothermal 
Field the level of dissolved minerals is rather low, and the increase in speci- 
fic gravity resulting from the solid phase can be taken as 0.01 - 0.02. In 
blowout prevention calculations, we general'y use the higher value (0.02). 
Accordingly, the specific gravity of the thermal fluid in the Yangbajing field 
is: 


ry, = 1.04 - 0.04 x 10 2t - 0.02 x 10 “t2 (4-2) 


5. The Stratigraphic Pressure Gradient as a Function of Depth in the Yangbajing 
Field 


Here we only discuss the pressure gradient laws for a liquid within the reser- 
voir. According to Figure 2, when the differential increase corresponding to 
any well depth is dZ, the increase in the stratigraphic pressure is: 


dP = O.ir, * dZ (5-1) 


















key! 4, etratigraphic pressure P 
b. depth 2 


Figure 2. Stratigraphic pressure calculation chart 


According to formula (}<l)a, we obtain dZ = 1/Kdt. We then substitute the values 
from equation (4-2), After integration and application of the following initial 
conditions: 

t = to and P = Po when Z = Zo, 


we obtain the following integral formula: 


10(?=P 0) ™ t + =? - t * tat + t? at = 
2 = Bo = 1.04 0.04 x “art x 10 0.02 x tht tot +t! y a9 (5-2) 


it can be shown that when the temperature to inside the well is between 120° and 
370°C, the calculation error resulting from replacement of the term t*+ ‘at + t’ 


in formula (5-2) by (‘sy is less than 0.5 percent. Then, letting 
tT. ta,* £ x 107? (5-3) 
we can simplify formula (5-2) to: 


10(P =~ Po) 2 . 2 7 
os 1.06 = 0.06T - 0.02 T ($-2)a 


the term JOC? = : ) om the left side of formula (5-2) is actually the specific 


gravity averaged over depth of the hot water in the well section between 7, and 7, 
and if we replace it by r 


te? 
. 0” -P ) 

“We ~ ee Z5 (5-6) 

then ty, = 1-04 - 0.00T « 0.02T? (5-5) 


in order to derive the stratigraphic pressure P at depth Z, we have only to 
convert equation (5-4): 


P= 0.1 ((10Po ~ r. 25) + £2). 


if we let: sH = 10 Po ty Zo (5-6) 


then P «0.1 (OH #4 Tua”? (5-7) 











In the above formulas, Po is the stratigraphic pressure (kg/cm’) within the 
reservoir at Che known depth Zo (meters); Zo is the depth of the known point 
(meters); P is the stratigraphic pressure at the calculation point Z% (ke/em*); 
% te the depth of the calculation pointe (metera); AHo is the pressure constant 
(in meters of a water column); ry, is the specific gravity averaged over depth 
of the hot water between dephts 79 and Z of the well section in the reservoir; 
and T is 1/100 the average temperature in the well section from Zp to 7, 


Below, we discuss the depth-averaged specific gravity Te and the pressure con- 
stant Ae of formula (5-7). 


A. The Depth-Averaged Specific Gravity 


We can see from equation (5<5) that the depth-averaged specific gravity between 
Zo and 7 is related only to the average temperature divided by 100 in that sec- 
tion. We know from the discussion of the function ¢(Z) that for most reservoirs 
(convective atrata), t is almost invariant with depth (with the exception only of 
transition zones which cold water is entering), and accordingly there is no neces- 
sary connection between the depth-averaged specific gravity and the depth. It is 
very important to recognize this fact, because it can greatly simplify the calcu- 
lation of the depth-averaged specific gravity. 


The above discussion was for single-stratum conditions. If one well passes 
through several reservoirs with different temperature gradients K or different 
temperatures, as in Figure 1, and if the known point 2, and the calculation 

point 7 are in different reservoirs, we rust first use equation (5-5) to find the 
depth-averaged specific gravities for the individual sections separately, b fore 
we can derive the average value. 


it should be pointed out that for ordinary cold water wells, the calculation is 
much simpler, since we can take "We equal to 1.0. 


5. the Pressure Constant 


The physical significance of the pressure constant can be seen from formula 

(5-6) and Figure 3. tt is equivalent to the height of the pressure head above 
the surface of the wellhead. This is a very useful constant which gives 4 vell's 
pressure-bearing state. 


The pressure constant is almost independent of well depth, and accordingly using 
it to describe the pressure-bearing state of a well is more universal than using 
the stratigraphic presgure or the pressure gradient, because separated from well 
depth conditions the latter two quantities are meaningless. 


10 














Figure }. Physical significance of the pressure conatant Alig. A tube filled 
with hot water of average specific gravity ry, and ite top sealed simulates a 
well with closed wellhead. The height to which the fluid in the manometer 
(eapecific gravity 1) rises (height of pressure head) is equivalent to the 
pressure constant Aly. 


Key: 

4. tube opening sealed; simulates well with closed wel ihead 
b. hot water (ry) 

c. U-tube manometer, water with \* l 

d. pressure Py 


The depth-averaged specific gravity and pressure constaits for certain welis in 
shallow reservoirs in the Yangbajing field are shown in Table 2. 


Table 2. Pressure Parameters of Some Welle in the Yangbajing Field 
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Methods of determining the location of the pressure constant: 


One of the methode given below can be used to determine the pressure constant 
AHoy the methods can be used for determinations under different circumstances. 
it te best to make determinations using several methods for mutual verification. 


(1) Actual Pressure Determination Method 


After a well which has been thoroughly blown out is closed, a downhole manometer 
is used at some depth 79 inside the reservoir to determine the pressure P, inside 
the wel , after which the calculation is made with equation (5-6). 


it must be pointed out that only pressure determinations made inside the reser- 
voir can be used for this calculation. To make the reading be ae close as 
possible to the real pressure, it ia best to lower the pressure sensor to a 
location where the fiuid passages are relatively good, for exampie a porous 
Stratum with rather large fissures or permeability. 


(2) balancing Mud Density Method 


When balancing occurs in the well drilling process, the specific gravity of the 
drilling mud can be finely adjusted so that the pressure of the mud coiumn and 
the stratigraphic pressure balance each other; the mud specific gravity ry, at 

the time of equilibrium and the bubbling or infiltration location Z» are recorded, 
and the pressure constant is calculated according to the following formula: 


AHy = (to = ty.) Zo (5-6)a 
(3) Wellhead Back Pressure Method 


When the location of a bubbling of infiltration stratum cannot be determined, 
this method may be used. This method too is used when the pressure coefficient 
SHo ie positive. It can be used to determine simultaneously the position 7, of 
bubbling or infiltration and the stratigraphic pressure Pp. 


Toe well ie pressurized with two kinds of muds with specific gravities r,; and r> 
which are both lower than that required to balance the stratigraphic pressure, 
and the casing pressures P-, and P.» for the closed well are determined sepa- 
rately, after which P» and Z»5 are calculated by the following formulas: 


since Po © Poy + 0.1 ry2_ 
and PrP, #® Pe? + 0.1 roZo 

~ (Pc; - Peo) 5-8 
accordingly, Zo oti (5-6) 


By subetituting Z_9 into any one of the formulas above, P, can be obtained. The 
value of Ay is then calculated by formula (5-6) using the values of Zo and Po. 








(4) dater Level Depth Method 





This method is used when there is 4 low pressure gradient, i.e., when the pres- 
sure constant ANo ie negative. 


Col¢ water ie introduced into the well, and the static water level 7; is deter- 
mined, after which the stratigraphic pressure Po can be calculated from the 
following formula: 


Po = O.1 ty (2p = 2) (5-9) 


where tr. = 1, the specific gravity of the cold water. The value of Alig can be 
calculated from equation (5-6) using Zo and Po. 


if the water infiltration level Z9 cannot be determined, two thin muds (specific 
gravity not greater than 1.10) with different specific gravities r, and ro, 
or plain water, can be introduced into the well, and the static level 2; and Zo 
determined, after which Z> cam be calculated from equation (5-10): 


since Po = 0.1 £) (Zo- 2) 
and Po = 0.1 rp (Zo- 29) 
» 12) - f222 
accordingly Zo r) - (5-10) 


it is easy to find Po from Zo, after which 4H9 can be found from equation 
(5-6). 


Below we derive formulas for the stratigraphic pressure gradient and the equiva- 
lent specific gravity. 


If we substitute the value of P from equation (5-7) into the pressure gradient 
formula and the equivalent specific gravity formula, we obtain: 


= - 
, Th Ty + a (t=-12) 


where GC is the geothersal pressure gradient and r, is the stratigraphic pres- 
sure equivalent specific gravity. Moreover, rm, = 106. G and rp have similar 
forms, and accordingly because of limited space we will discuss only rj. 
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Figure 4. Curves of rb(Z) for welis 6A and 13, Yangbajing field 


Key: 
a. (illegible) c. Well 13 
b. Well 6A d. depth, meters 


The curve for the stratigraphic pressure equivalent specific gravity for wells 
Nos 6A and 13 of the Yangbajing field, determined from equation (5-12), are 
shown in Figure 4. Well 6A is a well with a positive pressure constant, and 
well 13 has a negative pressure constant. The curves for the two wells have 
different directions of curvature. 


6. The Surface Patterns of the Stratigraphic "ressure Credient and the Zero- 
Pressure Surface Contour Map for the Yangbajing Field 


When choosing a well drilling process for new well sites and determining the 
surface distribution of the reservoir, it is necessary to understand the var '*- 
tion of the pressure gradient over the surface coordinates of the well loca- 
tions, i.e., the surface distribution law. 





Figure 5. Relationship of depth coordinate to elevation coordinate 
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For ease in comparing the parameters for different well locations, here we use 
the elevation above sea level inscead of the depth as previously. If we let the 
elevation be y at a location with depth Z, the wellhead elevation is h, as in 
Figure 5, we obtain: 


Zehe-y (6-1) 
and substituting into equation (5-7) we obtain 

P = 0.1 [Ao + rach - y)) (6-2) 
According to equation (6-2) there must exist some elevation above sea level hg 


at which the stratigraphic pressure P is equal to zero: y = ho, P = 0; substi- 
tuting into equation (6-2) we obtain 


ho = h + SMe 
Twa (6-3) 
and if we let dho = ho = h (6-4) 
we obtain Aho * OH (6-5) 
r 
Wa 


where ho is the elevation in meters at the zero-pressure point, and Aho is the 
difference in elevation between the zero-pressure point and the wellhead topo- 
graphy, and is called the excess altitude of the zero point. 


The sea level elevation ho of the zero-pressure point indicates the magnitude of 
the potential energy of the fluid; the excess altitude of the zero-pressure point 
Sho is similar to the pressure constant AHp in giving a well's pressure bearing 
characteristics, but AHo is in pressure units expressed as the height of a water 
column, while Aho is in geometric units. The value of hg is determined only by 
reservoir characteristics, while AHy and Aho are related not only to reservoir 
characteristics but also to the terrain elevation. With identical hoF, the 
greater the elevation of the wellhead, the smaller AHy and Sho, which completely 
explains why wells drilled in the northern part of the Yangbajing field are 
likely to infiltrate, while in the southern part they are likely to gush. 


By connecting a series of zero-pressure points we can form a zero-pressure sur- 
face. In order to express this zero-pressure surface on a flat diagram, we use 
the contour method, and call a map of this kind a zero-pressure contour map. 
Figure 6 is a zero-pressure contour map for shallow reservoirs in the Yangbajing 
field based on the zero-pressure point elevations of existing wells. 


The characteristics of some of the wells in the Yangbajing field are shown in 
Table 2. 


Below we given an example to explain the use of the zero-pressure contour map. 
If the surface position and wellhead altitude of well B which is about to be 


drilled is h = 4288 meters (see Figure 6), determine the specific gravity of the 
mud needed to balance the stratigraphic pressure in that well. 
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Figure 6. Zerc-pressure surface contours for shallow reservoirs on Yangbajing 
field 


Key: Large number beside well symbol is well number, numerator is zero- 
pressure point elevation, denominator is wellhead elevation. 


Answer: First we show the surface coordinates position of well B on the zero- 
pressure map (Figure 6), and read from the map the zero-pressure altitude ho, 

which is 4314 meters. Thus the excess elevation of the zero-pressure point is 
Sho = ho - h = 26 meters. 


On the basis of temperature data from well 5, we can use formula (5-5) to cal- 
culate the depth-averaged specific gravity ry . If the temperature is unknown, 


we can determine it approximately from the depth-averaged specific gravi- 
ties of wel's in the vicinity of well B (such as 6A and 9), obtaining Twa * 0.92. 


We obtain the pressure constant from equation (6-5): Aly = * Aho = 29.9 
meters water column. From equation (5-12), the relationship of the stratigraphic 
pressure equivalent specific gravity , and the depth for well B is: 


. OHy . 23.9 
r, ty, * 9 0.92 +5 
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In order to prevent gushing, the c -rect mud specific gravity t, should be 
greater than r, 0.02 - 0.04; if \ choose 0.02, then 


r = 0,94 + 2). 23.9 
m 


We can see from this equation that the specific gravity of the mud should be 
gradually decreased as the well becomes deeper. The maximum specific gravity 
is determined by the depth at which the top of the shallow reservoir is 
located; for example, if the depth of the reservoir Z is 30 meters, then the 
specific gravity of the mud used at the beginning of drilling should be: 

a9 oe 


.° 0.96 += = 1.74 


7. Use in the Drilling Process 


Here we will give a simple description of the use of the pressure gradient law 
in the problems described in section 1; a complete description is not the pur- 
pose of the present article. 


A. Estimation of the Pressure Bearing State of the Shallow Reservoirs in the 
Yangbajing Field 


A high-pressure constant indicates that the well is likely to gush, while in 
the reverse case it is likely to infiltrate. Accordingly, the potential for 
gushing (or infiltration) can be quantitatively estimated by means of the pres- 
sure constant AHpj. We can see from Table 2 and Figure 6 that the value of 

OH for the southeast portion of the field is higher than that for the north- 
west part, and on the mountain slope north of Zhongni Highway AHo even becomes 
negative. Accordingly, wells drilled in the southeastern section are likely to 
gush, while those drilled in the northwestern section are likely to infiltrate; 
in addition, we can also see that in the vicinity of a line connecting wells 2 
and 6A, the pressure constant is rather high, and the farther away to the side 
one moves from the line between, the lower AHg becomes. 


The area with high pressure constant values is in the vicinity of well 6A, and 
the maximum value is 24.1 meters (water column). On the basis of this datum, 

we are justified in beileving that the stable wellhead pressure value for wells 
drilled into the shallow reservoirs in the Yangbajing field that have been 
closed cannot exceed 2.5 kg/cm’. It is easy to calculate from equation (5-12) 
that for these wells (with a reservoir depth Z generally about 30 meters), we 
can control the stratigraphic pressure with a mud with a specific gravity of 
about 1.75 and prevent gushing. This completely answers the first question that 
was posed in section 1. Actual results have corroborated this conclusion. 


B. Selection of Well Suppression Method 
Geothermal wells have a rather high temperature, and the literature of many 


countries has reported the use of pressurized cold water cooling as a means 
of controlling gushing. This method was also used in the initial period of 
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drilling in Yangbajing. But under many conditions the method was ineffective. 
Analysis and calculations of the stratigraphic pressure gradient helped us to 
understand the essence of the phenomenon. We believe that the essence of cold 
water control depends on whether or not the pressure gradient of the cold water 
column can exceed the stratigraphic pressure gradient. This means that only 
when the stratigraphic pressure gradient equivalent specific gravity is less 
than the specific gravity of cold water can the pressurized cold water method 
achieve results. Otherwise, it is necessary to use the mud well suppression 
method or stay with the drill rod suppression method cc the capping suppression 
method. 


Here we must point out that gushing and infiltration are two expressions of the 
same essence: the former is violent and the latter gradual. In general when 
the wellhead excess pressure is maintained below 0.5 kg/cm* it is the latter 
form that occurs. For some wells with rather low pressure constants (but 
larger than zero), if a sufficient amount of continuous cold water injection 

is used to maintain the liquid column pressure steady, a relatively low well- 
head excess pressure can be maintained in a constant state, at which time the 
goal of prolonged bubbling but no gushing can be achieved. But this "well 
suppression" is incomplete. 


C. Selection of the Length of the Uncased Bottom Hole and the Position of the 
Casing Point 


We can see from formula 5-12 that when the pressure constant AHp is positive, 
the stratigraphic pressure equivalent specific gravity tr, decreases hyperbol- 
ically as the depth Z increases. The shallower the well, the more rapidly it 
drops off. Accordingly, when the length of the uncased bottom hole is too great, 
the difference in the mud specific gravities required at the two ends of the un- 
cased section is increased, with a high specific gravity required at the top and 
a low specific gravity required at the bottom, so that it is difficult to use a 
mud with one specific gravity to satisfy the two different requirements. This 
is the reason why geothermal wells are extremely likely to show upward gushing 
above and downward infiltration below. This phenomenon is even more pronounced 
when drilling into shallow strata. The method of overcoming it is to select 

the location to which the casing penetrates in accordance with the infiltration 
suppression pressure gradient and the stratigraphic pressure gradient, to control 
the length of the uncased bottom hole, and at the same time to use an infiltra- 
tion plugging method. In addition, the mud specific gravity should be correctly 
calculated for different depths. 


D. Selection of a Method of Inducing Gushing 


The pressure gradient data are straightforward indications of the degree of 
difficulty of inducing gushing of a well after it has been completed. Wells 
with a pressure gradient greater than 0.1 can conveniently be induced to gush 
by decreasing the mud specific gravity; those with a gradient less than 0.1 can 
only be induced to gush by gas lift or pumping of water. Some wells in which 
the pressure gradient is close to that of plain water can be induced to gush by 
the “heating effect," which consists of introducing plain water and waiting for 
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it to become heated, when its specific gravity will decrease, thus inducing 
gushing. The stratigraphic pressure gradient data are very useful in select- 
ing and calculating the above flow-inducing methods. 


In other countries the heating effect method after introduction of cold water 
is extensively used for high temperature geothermal wells with temperatures 
between 200° and 300°. But it must be seriously pointed ou. that this method 
is generally useless for certain negative-pressure-constant wells with a tempe- 
rature below 200° and with cold water runoff in the shallowly buried section. 
Under these latter circumstances, the saturation pressure of the hot fluid is 
rather low, and there is no way of overcoming the pressure of the water column 
(with a rather low temperature) in the well. 


E. Selecting the Drilling Fluid 


In order to control gushing, prevent infiltration, and also to prevent the 
reservoir from becoming clogged, the best approach is to use the pressure 
balancing method of drilling. The stratigraphic pressure gradient data are 
useful in selecting the kind of drilling fluid to be used. In general, when 

the stratigraphic pressure equivalent specific gravity is equal to or gieater 
than the specific gravity of plain water, plain water or a drilling iluid of the 
appropriate specific gravity can be used to drill the well, but when the equiva- 
lent specific gravity is less than that of plain water, foam, fog or air can be 
selected for drilling, depending on the magnitude of the equivalent specific 
gravity. The low specific gravity fluid method of well drilling is promising 
for wells in the northwestern part of the Yangbajing field and in wells whose 
deeper sections have low pressure gradients. 


8. Conclusions 


A. The stratigraphic classification used for oil or gas wells should not be 
used to classify strata of geothermal fields. This article proposes an approach 
to pressure classification. 


B. The relationship between the stratigraphic pressure gradient of a geothermal 
field and depth can be described by a hyperbolic function. In the shallow sec- 
tion, the parabolic curve changes steeply, while in the deep section it slopes 
more gently. Accordingly particular care must be taken when drilling into 
shallow reservoirs. 


C. The direction of curvature for wells with positive and negative pressure co- 
efficients are opposite, but they all have the depth-averaged specific gravity 
of the heat fluid as their limit. 


D. A method of drawing a zero-pressure surface contour map is proposed. It is 
helpful in estimating the stratigraphic pressure gradient situation at different 
surface positions. 
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Afterword 


Although this article is intended to apply to high-temperature moist steam 
geothermal fields, its principles are also applicable to hot spring wells, cold 
water wells, exploratory drilling in solid mines and other engineering drilling 
exploration. Under these conditions, because the temperature is rather low, 
(for the most part normal) the depth-averaged specific gravity of water Ty, can 
be taken as 1.0, so that the various calculations become much simpler. 
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APPLIED SCIENCES 


HAJIF-1 AIRCRAFT STRUCTURAL ANALYSIS SYSTEM DESCRIBED 
Beijing GUO.1 HANGKONG [INTERNATIONAL AVIATION] in Chinese No 2, 5 Feb 80 pp 2-4 


[Article by Aircraft Structural Analysis System Development Group: "Brief Intro- 
duction to the HAJIF-1 Aircraft Structural Analysis System") 


[Text] The HAJIF-1 (Hangkong Jiegou Fenxi Xitong 1) Aircraft Structural 
Analysis System is a large-scale general-use set of finite element structural 
analysis computer programs independently developed by our country. A study 
lasting 10 months has shown the system's operation to be satisfactory, so that 
it is now ready for formal use. 


Back in the late 1950's, our country began to use computers for stress and defor- 
mation analysis of aircraft structures. At the time, a displacement wethod using 
beams and torsion resistance boxes was used. In the early 1960's a standardized 
program was developed using a direct stiffness method and substructure techniques. 
Starting in the 1970's, the limited-element method was used to develop many 
specialized and general programs for use in structural analysis and summariza- 
tion. In 1976, the China Aviation Research Institute (CAE) decided to bring 
together groups belonging to four of its research institutes in a special devel- 
opment group in order to begin development of a large-scale general structural 
analysis system; it started by developing the HAJIF-1 linear static analysis 
section. 


Design Approach 


In order to make the new system easier to use, more efficient and broader in 
application than previous structural analysis programs, the development work 
concentrated on automatic generation of starting data, a special user-oriented 
structural analysis language, and a highly efficient, high-order sparse matrix 
analysis method. The specifics of the design approach include: 


l. Generality of application. A multilevel substructure analysis technique was 
used, making it possible to analyze any type of large-scale complex structure. 


2. Flexibility. We used a user-oriented structural analysis language which 


could carry out 15 fixed routines contained in the system and could also util- 
ize specially-defined operations written by the user. 
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3. High efficiency. The system features specialized data statements and 
automatic generation of starting data, and offers the user two alternative 
solution methods. 


4, UExpandability and ease of revision. The system has a modular structure, 
so that to revise or add new capabilities it is only necessary to revise or add 
the relevant module, without altering any other modules. 


5. Interrupt and diagnostic features. The system has an interrupt processing 
capability, a step-through capability, and can restart after a present inter- 
rupt is found. It has a certain error diagnosis capability, and it can emit a 
warning signal or stop operation depending on the nature of the error. 


6. Easy machine jnterchangeability. With the exception of a small part of the 
program which is written in assembly language, the program is written in FORTRAN 
IV, which makes it easy to use it on computers of different types. 


System Capabilities 


The HAJIF-1 system can carry out linear deformation and stress analysis of thin- 
walled structures, crossbeam frames and framed girders under external load, and 
will soon have the additional capability of analyzing three-dimensional continu- 
ous elastic bodies. There are no limitations on the types of problems that can 
be solved. But because the computer on which it is currently being used has a 
relatively small memory (256 Kbyte) and is rather slow (fewer than 100,000 ope- 
rations per second in double-precision floating point arithmetic), the scale of 
the problems is limited to a maximum of 10 substructure levels, and a maximum of 
99 individual substructures, with no more than 3,000 degrees of freedom for each 
substructure and no more than 6 degrees of freedom for each node. The system 
has 30,000 FORTRAN and assembly language statements. 


The system offers 15 fixed routines, including triangularization analysis, dis- 
placement, stress analysis and influence number calculation based on assembly of 
a stiffnese matrix; triangularization analysis, displacement, stress analysis 

and influence number calculation based on assembly of a boundary stiffness matrix; 
displacement; stress calculations; influence number computation, and the like. 


The system's elements store contains 18 elements including various rod elements, 
beam elements, plate elements, thin surface elements and rigid elements. The 
materials store contains 7 materials, including LY12 (plate and shaped stock), 
LC4 (plate and shaped stock), 30CrMnSiA, 30CrMnSiNi2A and 1Crl8Ni9Ti. 


Mechanics Design of the System 


The HAJIF-1 system divides the nodal degrees of freedom into five types: degrees 
of freedom of nodes within substructures; degrees of freedom of nodes at sub- 
stivucture boundaries; degrees of freedom at moving boundaries, primarily used 
when the generalized displacement in the direction of one degree of freedom is 

a fixed value; permanently fixed degrees of freedom, primarily providing con- 
straints when eliminating motions of structural rigid bodies; and symmetrical 
and antisymmetrical degrees of freedom, primarily used in symmetricity treat- 
ment of structures. 


22 








The HAJIF=<1 system provides five rectangular coordinate systems: substructure 
courdinate systems; element coordinate systems; nodal coordinates, primarily used 
in eliminating abnormal states arising from the formulas as a result of very low 
strength in the direction of one of the degrees of freedom; assembly coordinates, 
principally used in providing parameters for standard-geometry assemblies of etruc- 
tures and automatic generaiion of the calculation model grid; and structure co- 
ordinates, which function as common reference numbers for all of the structures. 


The basic steps of the solution by HAJIF<-1 are: assembly of the eubstructure 
stiffness matrix, handling of constraint elements, substructure analysis, symme- 
try processing, interconnection of substructures, and multilevel substructure 
anaiyels, displacement, stress computation, and influence number computation. 


Syetem Mathematical Methods 


In the static analysis. trisngularization analysis of the stiffness matrix 
occupies the largest part of computation time. Accordingly the efficiency of 
the high-level sparse matrix solution sethod very strongly affects system effi- 
ciency. HAJIF-1 offers two solution methods which can handle sparse matrices 
of any widt) and are very efficient. 





1. Automatic partitioning and effective column solution method, which uses two 
techniques to make it efficient. It avoids all empty operations which do not 
change the calculation results, and minimizes intermediate multiplication and 
addition results in the inner cycles; and depending on the size of the work area 
in the computer internal memory, it divides the columns of the matrix it contains 
into two types: those in the partitions, which are the principal columns, and 
those included among the effective columns, which are those whose influence: 

must be calculated. During each step, only the elements needed for that step are 
entered into memory, which decreases the number of transfers between internal and 
external memory. 


2. The superelement matrix analysis method, which simultaneously divides the 
stiffness matrix into many original level submatrices, and arranges these sub- 
matrices on-the basis of their corresponding positions in the stiffness matrix 
#0 a8 to achieve a two-level direct-indexed matrix. If the order of the second- 
level address-indexed matrix is still too high, it can be further divided using 
the above method to obtain third-level addressed indices. This data organiza- 
tion method is relatively well suited to third-generation computers, and makes 
tranefer of data in ‘nternal and external memory convenient and increases effi- 
ciency. 


Program Organization 


The HAJLF-1 program organization is on the two-level principle (Figure 1). The 
first level ia the system individual capability modules, which are independent 

of each other, and merely carry out their individually specified functions, with- 
out being able to affect each other. The second level is the supervisor program, 
which is used to control the various functional modules, and to control the en- 
tire calculation and analysis process. The supervisor program is resident in 
main memory, while the various functional models are rolled in and out. 
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Figure 1. Overall structure of HAJIF-1 





Key: 
JKCX: Supervisor program SJSH: Data generation module 
JSMK: Interpreter module YXSJ: Initial shape data 
YONG HU...: User manual THLX: Geometrical snapes store 
WSCL: External processing YSWZ: Element position store 
JPRP: Assemble stiffness matrix DPY: Equivalent load displacement 
YLJS: Stress computation store 
YSKU: Element store DPJS: Equivalent load computation 
FIKP: Analyze stiffness matrix PDSH: Unit load formation 

and loads YXJS: Influence number computation 
ZJZP: Substructure assembly DY: Printout module 
WYHD: Displacement SCYL: Print stress curve 


[he system provides a user-oriented structural analysis Language in order to 
adapt the multiplicity of structures and the different requirements of differ- 
ent usere. This Language consiste of 92 statements which use Chinese romani- 
ration letters. The system interpreter module translates the initial program 
weitten by the user in the structural analysis language into a string of inter- 
nal code which can be recognized by the system supervisor program. The super- 
visor program then organizes the computation process on the basis of the 
information contained in the internal code. In addition to the 15 above- 
mentioned routines provided by the system, the user also may use the structural 
analysis language to write his own subroutines, thereby meeting the need for 
tlexibility. 


Data Generation in the System 


HWAJIF-1 hae a data generation module which can automatically generate the start- 
ing data, greatly simplifying the user's data preparation work. Data generation 
can be civided into external load calculation and estructural data generation. 


The task of external load computation is to assign the distributed aerodynamic 
forces and inertias equivalently over the structural mechanics model's nodal 
points and assure that the entire aircraft's external serodynamic loads are in 
comp ete equilibrium with ite inertial loads. WAJIF-1 can calculate the aero- 
dynoemic force distribution over the wings, tail and body on the basis of wind 
tunne. results or strength standards, can calculate the inertial force distribu- 
tion, and of course cau. use loads specified by the operator. 
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The structural data generation module can, on the basis of user-specified con- 
trol pointe and section lengthe, automatically carry out a grid analysis of the 
standard geometric components: lines, curves, planes, curved surfaces, and 
form the corresponding nodal coordinates. The data generation module provides 
a total of 13 different geometric assembly shayes, including standard curvi- 
Linear segments, composite curvilinear segments, line segments, arc segments, 
parabolic segments, conic sections, parabolic sections, four-sided planar sec- 
tione bounded by 4 composite curves, irregular nodal geometry component types 
and the like. 


For element node information, the data generation model provides five element 
position generation statements. 


For conventence in processing large quantities of data that continually reappear 
according to a specific pattem, the data generation model contains two constant 
difference data repetition statements and one data repetition statement. In 
summary, it can be seen that the various steps in HAJIF-1 from load calculation 
and structural data generation to calculation of structural deformation and 
stress distribution are connected into a whole, and accordingly they are con- 
venient for the user and efficient. 


An Example of System Utilization 


Atter overall system debugging was completed, a total of 26 sample problem cal- 
culations were made. These problems ranged from sma'l scale to large scale and 
from simple to complex, and finally a complete stress analysis of a certain air- 
craft was carried out. The entire body of the aircraft was divided into 20 sub- 
structures ag shown in Figure 2. The largest substructure had 2,238 degrees of 
freedom, and the whole aircraft, with 27,446 elements, had 30,502 degrees of 
freedom. There were 9,188 input data cards. The multilevel substructure tree 
is shown in Figure 3}. 











Figure 2. Substructures of a certain aircraft 
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Figure 3. Multilevel substructur® tree 


Kev: 

a. Fitth level substructures a Second levei substructures 
b. Fourth level substructures e. First Level substructures 
ec. Third level substructures f. Substructures 


CWKZ (1) Indicates performance of error diagnostics 
SJSH Data generation 


KOSX 1/2, 3/4, 5/6, 7/8, 9/10, 11/12, 13/14, 15/16, 17/18, 
101 (1, 2), 102 (3, 4), 103 (5, 6), 106 (7, 8), 105 (9, 10), 
196 (il, 12), 107 (13, 14), 108 (15, 16, 17, 18), 109 (19, 
20), 201 (101, 102, 103), 202 (105, 106, 107), 301 (202, 109), 
401 (201, Wl), 501 (401, 108). Tree structure for the calcula- 
tion; 1/2 indicates that substructures | and 2 are symmetrical, 
and 102 (1, 2) indicates that first-level substructure 0] con- 
sists of substructures 1 and 2. The rest is analogous. 





1.05 Fixed routine proceeding from assembly of stiffness matrix 
to stress computation. 


DY (ZOLYSYL, 100; ZO2YSYL, 101, ZO3YSYL, 102, ZO4YSYL, 103) 


Printout: stress on first substructure; elenent stress for 
second substructure with ree‘dual strength coefficient 

n >0.95; element stress for third substructure where »> 1.3; 
stress, 1, and nodal force on components of fourth substruc- 
ture with maximum residual strength. 


AAN Program completed. 


HAJIP-l's capabilities are limited to solving linear static problems at constant 
temperature. Such problems as dynamic analysis, stability analysis and linear 
influence analysis of materials and geometry will be solved in later versions of 
the structural analysis system. HAJIF-1 will be transferred to a larger and more 
modern computer. 
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Details of the aircraft estructural analysis system can be found in the "Deacrip- 
tion of the Aircraft Structural Analysis Syatem, Version 1 (HAJIF-1)," in the 
the “Aircraft Structural Analysis System Version 1 (HAJIF<1) User's Manual,” and 
in a succession of articles which will appear in aviation journals and conferences 
on computational mechanics. 
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TITLe+: “Supersonic Seil Wing" 


SOURCE: Beijing HANGKONG ZHISHI [AERONAUTICAL KNOWLEDGE) in 
Chinese No 11, Nov 80 p 14 


ABSTRACT: This article introduces a new wing design which can 
provide improved performance of a supersonic aircraft by utilizing 
favorable shock wave interference effects. In particular, the 
physice of shock wave formation over a wing during supersonic 
flight is explained. A specific sail wing design developed by 

the MeDonnell Douglas Company is introduced and a tutorial 
discussion is presented to illustrate how the lift-to-drag 

ratio of the wing can be enhanced by taking advantage of the 

nigh pressure region below the wing surface which is created by 

an oblique shock wave emanating from the aircraft nose. 


AUTHOR: YUE Min (1471 2404) 
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TITLE: “Infrared Homing Guidance” 


 ~OURCE: Beijing HANCKONG ZHISHI [AERONAUTICAL KNOWLEDGE) in 
Sninese No 11, Nov 80 pp 20-21 


ABSTRACT: Infrared homing guidance is a passive guidance tech- 
nique most commonly used on air-to-air missiles. It uses optical 
sensors which can seek targets that emit infrared radiation. In 
general, an infrared homing device employs the proportional navi- 
gation guidance scheme which steers the missile according to the 
measured angular rate of the line-of-sight rotation between the 
target and the missile. An infrared homing system has the follow- 
ine components: 1) optical system; 2) photo-sensitive detector; 

') signal voltage modulator; 4) electronic circuits; and 5) track- 
ing mechanism. Current research efforts in infrared homing gui- 
dance are concentrated in the development of highly sensitive 
detectors, fast-response tracking mechanisms, and automatic search 
and track capability. 
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TITLE: “High Altitude Simulation Test Platform for Aircraft 
Turbo Engines” 


SOURCEs Beijing HANGKONC ZHISHI [AERONAUTICAL KNOWLEDGES) in 
Chinese No 11, Nov 80 pp 28-29 


ABSTRACT: A high altitude simulation test platform is designed 

to create the pressure and temperature conditions of high altitude 
environment in order to test the performance of aircraft engines. 
Specifically, its functions are to: 1) certify and modify engine 
performance; 2) provide performance curves for airplane design; 

3) eliminate ongnee troubles during flight; and 4) provide new 

data for high altitude engine research. In this article, the basic 
units of a high altitude simulation test platform are introduced 
and the simulation process of creating the stagnation temperature 
and pressure conditions of engine inlet and the static pressure 
condition of engine exhaust are explained. The key components of 
a test chamber which include the front chamber, the flow regulating 
screen, the inlet nozzle, the rear chamber, and the thrust measur- 
ing device are described. 


AUTHOR: WANG Chen [3769 2525] 
ORG: None 
TITLE: “MATRA Super 530 Air-to Air Missile” 


SOURCE: Beijing HANGKONG ZHISHI [AERONAUTICAL KNOWLEDGE) in 
Chinese No 11, Nov 80 pp 32-34 


ABOTRACT: The MATRA Super 530 is a new French air-to-air missile 
to be equipped on the MIRAGE F-1 fighter airplane. It is designed 
specifically to intercept high altitude, high speed airplanes such 
as the MIG-25. The MATRA Super 530 was developed on the basis of 
ite predecessor, the MATRA R530 missile and was first deployed by 
the French Air Porce in December, 1979. The operation of a MATRA 
“uper 530 miesile can be either manually controlled with the aid of 
a fire control radar or automatically controlled by an on-board 
computer. The main advantages of the MATRA Super 530 missile are 
as follows:1) good aerodynamic characteristics with high degree of 
mobility, stability, and high speed (Mach 4-5) capability; 2) highly 
sensitive idance system with high-gain antenna, guidance computer 
and auto-pilot; 3) solid propellant engine to provide high thrust- 
to-weight ratio; 4) warhead with a kill radius of over 20 meters 
and a radio-controlled fusing device; 5) light-weight honey comb 
steel structure which can withstand high temperature created by 
aerodynamic heating. 
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AUTHORS; CHEN Tao [7115 3447] 
WU Yenyu [0702 3601 3768) 
ZHANG Xuexun [1728 1331 6104) 
KUNG Olnaxin [1313 1987 2450) 


ORG: All of Institute of Forestry and Pedology, Chinese Academy of 
Sciences 


TITLE: “Amelioration of Cadmium Soil in Chanqshi Irrigation Area and 
Prevention of Cadmium Contamination in Rice" 


SOURCE: Beitding HUANJING KEXUE [ENVIRONMENTAL SCIENCE] Vol 1 No 5, 
10 Oet 80 pp 7-1) 


ABSTRACT; The study on how to ameliorate the cadmium soil in the 

Zhanashi irrigation area of Shenyang, and prevent cadmium contamination 

in rice concludes that: 1) lime is the best amendment in terms of 

its neutralization of acidic soils, inhibition of rice absorption of 
cadmium from soils irriqated with acidic heavy metallic sewage, and 
economy; 2) ph is the main factor that affects rice absorption of cadmium, 


[Continuation of HUANJING KEXUE Vol 1 No 5, 30 Oct 80, pp 7-11) 


which drops against rising pH; and 3) rice absorption of cadmium is mainly 
at the later stage of growth, and therefore, more frequent irrigation at 
the earing stage is recommended for the prevention, 


The authors express thanks to Vice-director Gao Zhengmin [7559 2163 3046] 
of the Institute, and Professor Yao Guigeng [1202 2981 5087] of Shenyana 
Agricultural College for reading and correcting the manuscript. 





AUTHOR; TIAN Bingshen [3944 3521 3947) 


ORG; Atmospheric Chemistry Section, Gansu Institute of Environmental 
Protection 


TITLE; “Preliminary Study on Photochemical Smog Pollution at the Siqu 
District of Lanzhou" 


SOURCE: Beijing HUANJING KEXUE [ENV!®ONMENTAL SCIENCE) Vol 1 No 5, 
10 Oct 80 pp 24-30 


ABSTRACT; The study on the photochemical smog pollution over the Siqu 
district of Lanzhou shows that the frequency of such photochemical smog 
derived from nearby heavy and chemical industries can be higher than 
35%, and the degree of pollution is indeed serious. The causes of 

such photochemical smoq pollution is attributed to the strong solar 
radiation through the atmosphere th*t contains definite concentration 
of nitrogen oxides, crabohydrides and sulfur dioxide, Contrary to 
those at Los Angeles, the smog here is caused mainly by discharged 
waste gases from nearby industries, 


The author expresses thanks to the Department of Environmental Chemistry 
of the Beijing University for direction and assistance, and to Prof, 
Tang Xiaoyen [0781 1321 3508) for reading tue manuscript. 


AUTHOR; CHU Jiachenag [5969 1367 2052] 


ORG: Water Transportation Department of Ministry of Communication's 
Academy of Sciences 


TITLE: “Preparation of a Dispenser for Oil Spills” 


SOURCE: Beijing HUANJING KEXUE [ENVIRONMENTAL SCTENCE] Vol 1 No 5, 
10 Oct 80 pp 62-66, 47 


ABSTRACT: Based on a study on the characteristics of HLB values and 
surface active aqents in the compounding of highly efficient and low 
toxic oil] spill dispensers, the author reports the preparation of two 
such preparations, with the adoption of such low toxic emulsifying 
agent as the ester type non-ionic surface active agent and paraffin 
hydrocarbon solvents. The HLB values for best emulsifying efficiency 
was found to be a tertiary compounding of £0(4), EBO(11) and FO(6) 
Oleats, but subsect to the surface active agents composed of 
polyethylene oleate and polyethylene sorbitolate at appropriate 
proportions, The properties of the preparations with regard to their 
emulsifying power, biological toxicity, biological breaking, flash 
point and viscosity are finally evaluated, 








AUTHORS; CHEN Guohua [7115 0948 5478) 
YANG Zongyu (2799 1350 3022) 


ORG; Both of Southwestern Institute of Applied Magnetism, Chinese 
Academy of Sciences 


TITLE; “Application . High Gradient Maqnetic Separation to Water 
Pollution Control" 


SOURCE; Beijing HUANJING KEXUR [ENVIRONMENTAL SCTENCE] Vol 1 No 5, 
30 Oct 80 pp 71-75 


ABSTRACT; Advocating the advantages of the modern application of hidh 
gradient maqnetic separation to water pollution control, the authors 
introduce the basic principles of the physical phenomenon of high 
gradient maqnetic separation, the structure and separation of high 
qradient maqnetic separators, and its application to the treatment of 
sewage from iron and steel industries, urban water works, and cleaning 
of drinking water, It is anticipated that adoption of superconductor 
technique will further increase the magnetic gradient, decrease power 
consumption, extend the separation range, and reduce cost, that water 
circulation will be further perfected, while river and lake pollution 
eliminated, 
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AUTHOR: LIN Chao (2651 6389] 
LI Changwen [2621 2490 2429] 


ORG: Both of the Department of Geography, Beijing University 
TITLE: "Studies on Land Types in Mountain Regions North of Bei jing" 


SOURCE: Beijing DILI XUEBAO [ACTA GEOGRAPHICA SINICA] in Chinese Vol 35 No 3, 
Sep 80 pp 187-199 


TEXT OF ENGLISH ABSTRACT: The mountain regions north of Beijing occupy the western 
part of Yan Shan which runs eastward to Shanhaiguan and forms the northern border 
of the North China Plain. The climate belongs to the semi-humid warm temperate 
type, with a cold winter and a hot summer. Average annual precipitation is about 
800 mm in the front ranges and decreases to about 500 mm in the sheltered valleys. 
Deciduous forest of a semi-xerophytic type is the common vegetation. Drap soil is 
developed under this climatic and botanical environment. 


Field investigations revealed that a great variety of landscape exists in this 
mountain region, A system of classification of land types has been attempted. 
This classification is based on the Knowledge of physical geography and especially 
on the analysis ot the physical factors controlling the differentiation of the 
landscape. 


First, altitudinal zonation has to be considered. Although the altitude is not 


[Continuation of DILI XUEBAO Vol 35 No 3, Sep 80 pp 187-199] 


very high, from about 50 m in the foothills to over 2,000 m at the summit, vertical 
zonation of natural landscape is clearly manifested, Three vegetation-soil zones 
may be distinguished: 

l. The lower part of the mountain, below 1,200 m, is covered by shrubs and grasses 
of semi-xerophytic and meso-xerophytic types. Small patches of Pinus tabulaeformis, 
Quercus dentala and Fraxinus bungeana may be found here and there in sheltered and 
protected areas. This is a secondary growth, the result of interference of human 
activities since historic times. Deprived of the protection of forest cover, ero- 
sion is very intensive on this part of the mountain. Mass wasting, such as mud- 
flows, isa common feature which causes great damage to property and loss of life. 
2. Above 1,200 m and up to 1,850 m is the second zone. It is rather well covered 
by a deciduous forest, dominated by various species of oak (Quercus aliena, Q. vari- 


abilis, Q. accussima), birch (Betula piatyphy lia, B. duhurica) and poplar (Popular 
davidiana). Shrubs and grasses form the undergrowth of the deciduous forest, but 


they are different species from those in the lower mountains. Brown forest soil is 
the typical soil. Due to the protection of the forest cover, mudflows seldom occur 
in this zone. 

3, Above 1,850 m to the top of the mountain the deciduous forest is superseded by 
mountain meadow which is underlaid by a layer of black meadow soil. It constitutes 
a distinctive zone by itself. The meadow is used for pasture during the summer. 








Next to altitudinal zonation, lithology deserves special attention. Lithological 
character is an important factor in the formation of landforms, soils and 
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vegetation. For example, we found granite and limestone forming a very different 
landscape in this part of Yan Shan, Granite usually forms rounded or undulating 
landforms, covered with a thick layer of regolith and a mixed forest of pine and 
oak, Brown forest soil is developed wherever granite occurs. In Limestone areas 
karstic Landforms are found, with only a thin veneer of weather materials. Cypress 
(biota orientalis) usually takes the place of pines. Drap soil is usually developed 
wherever limestone occurs. Loess which appears in different altitudes forms a 
typical landscape by itselt, characterized by flat or undulating top and vertical 
slope under water erosion, The typical soil is drap soil, most of the loess- 
covered land being cleared for cultivation. Alluvial deposits along the valley con- 
stitute another type of landscape. The new deposits beside the river channel form 
flocd plains, while the older deposits further away from the river form terraces. 
The flood plains and terraces are the most fertile land in the mountain region and 
become the center tor agriculture and settlement. However, the occurrence of summer 
tloods may be a hazard to the crops and inhabitants in the valley. 





Mountain Landscape is further diversified by the erosion of running water. Mountain 
slopes are intensively dissected and form mesorelief and microrelief. Slopes may be 
gentle or steep, convex or concave. The effect of aspects of slopes is also quite 
evident. In the upper part of the deciduous forest zone, south facing slopes are 
usually covered by shrubs and grasses and only the north facing slopes are forested 
with deciduous trees. In the lower mountain zone, the south facing slopes are 
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usually cleared for growing fruit trees while the north facing slopes are forested 
with pines. 


based upon the analysis mentioned above, the land units are organized into a hier- 
archy of descending, size as tollows: 

Order I. (Land system or mestnostch), altitudinal zone with common type of climate, 
vegelation and soil, without regard to geological and geomorphic features. Litho- 
logical and geomorphic differences may form a sub-order within the zone. 

Order Il. (Facet or urochishche), a part of the altitudinal zone formed by grouping 
of different facies, having common mesorelief, rock materials, drainage, vegetation 
and soil, 

Order [II. (Facies or element), simplest unit of landscape, with homogeneous 
lithology, microrelief, drainage, microclimate, vegetation and soil. 


All the units are recurrent in their distribution within the respective order. The 
individual units (local forms) are grouped into abstract forms (species) according 
(o their homogeneity or similarity. Species of land units may be further grouped 
into a-scociation. 


the procedure of mapping and classification of land units may proceed from lower 


order (sm iller size) to higher order by grouping of the smaller units into larger 
units or, vice versa, by subdividing the larger units into smaller units. 
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Mapping of land types \ carried out on different scales. For the larger units, 
®.@+, the land systems, . map on the scale of 1:200,000 has been produced, This is 
suitable for the use © he provincial authorities. For the facets, a map on the 
seale ot 1:100,000 ha | ‘en produced. This is intended for the use of the county, 
A few waps for distrit ion to land units of the lower order (facies) have also 
been prepared on the scale of 1:10,000,. These may be used by the commune. It has 
been found that not all facies can be shown even on the scale of 1:10,000, so some 
facies have to be left out or combined. 


in this report the abstract land units are described, with descriptions on area, 
distribution, landforms, water conditions, soil, vegetation and land utilization, 
which would be convenient for reference, 


AUTHOR: ZHONG Gongfu [6945 0501 3940] 
ORG: Guangzhou Institute of Geography 


TITLE: "Mulberry-Dyke--Fish-Pond on the Zhujiang Delta--A Complete Artificial 
Ecosystem of Land-Water Interaction" 


/OURCE: Beijing DILI XUEBAO [ACTA GEOGRAPHICA SINICA] in Chinese Vol 35 No 3, 
Sep 80 pp 200-209 


(EXT OF ENGLISH ABSTRACT: The “Mulberry-Dyke--Fish-Pond" system is a peculiar farm- 
ing method developed by the working people as a result of making full use of land 
and water resources in their longtime farming practice in the Zhujiang Delta. It is 
an interrelated system which brings into tull play the production potential and 
wakes reasonable use of land-water resources. It promotes the development of dif- 
terent branches of agriculture and its related processing industry. "“Mulberry-Dyke-- 
Fish-Pond" is a complete scientific ecosystem, 


Mulberry -Dyke--Fish-Pond is characterized by its complex and diversified cycle. 
Mulberry growing, sericulture and fish breeding are interrelated. They promote each 
other and their circulating operation is complex and diversified. In terms of the 
tood chain, the whole operation follows the following pattern: mulberry leaves are 
fed to silkworms, silkworm excreta is fed to fish, and mud from the fishpond is used 
as good manure for mulberry trees. This circulating process begins with mulberry 








[Continuation of DIL] KUEBAO Vol 35 No 3, Sep 80 pp 200-209] 


growing, then moves on to sericulture and ends up with fish farming. In this closely 
related process, each link has either a positive or adverse effect on the others. 
There is a popular saying in Shunde which goes something like: The more Luxuriant 
the mulberry trees, the stronger the silkworms, the fatter the fish, the richer the 
pond, the more fertile the dyke and the more the cocoons. This illustrates in full 
the circulation relation of Mulberry-DykKe--Fish-Pond. In this system the mulberry 
tree is the producer making use of sunlight energy to produce mulberry leaves, the 
silkworm is the first consumer eating mulberry leaves, the fish is the secondary 
consumer having silkworm excreta (and pupas) as food, and the plankton in the pond 
is the reduction agent bringing the organic substances, such as N, P and Ky, back to 
the pond. 


The producing cycle of mulberry growing, sericulture and fish farming has as its 
base (he circulation of substances in soil, water and living beings between the pond 
and dyke. The relations of mulberry trees, silkworrs and fish to dyke, pond and 
atmosphere are those of substance-energy exchange and transformation between living 
beings and environment. Any change in morphology, productivity, quality, etc., dur- 
ing the growing process of mulberry trees, silkworms and fish is the result of 
interaction between mulberry trees, silkworms, fish and elements of the environment, 
such as temperature, light, water, air, soil, etc. 


With an interdependent relationship between the living beings and environment, the 
Mulberry -Dyke--Fish-Pond is a complete ecosystem in which substance circulates and 
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energy is exchanged and transtormed. It is a scientific farming method as well. In 
various cycles, big or small, the mulberry is the base and the fishpond is the key 
link. The pond can feed fish as well as the dyke. With the pond being taken good 
eare of, there can be high yield in fish, and that also promotes the output of 

other crops. 


‘he role of the Mulberry-Dyke--Fish-Pond is: 

(1) Bringing into full play the potential of the local resources. Located in the 
southern subtropical zone, each vear large quantities of mulberry leaves are col- 
lected, enough to teed eight to nine batches of silkworms, hence eight to nine batch- 
es of cocoons. Fish fry from the middle and lower parts of the Xijiang River are 
wade use of in the development of pond fish farming. The river water and tidal 

water are made use of in self-control drainage and irrigation all year round. Mud 
‘rom rivers and ponds is used as an important manure. Silkworm excreta is used to 
feed tish. Thus the land is being made full use of, 


(2) Abundant sources of income and huge economic gains. The average income from 
every mu of mulberry trees (including that from mulberry leaves, cocoons, fish, 
veget les, dry land crops and mulberry branches) is 400-500 yuan. If cocoons are 
exported, income will be increased, about four times that of rice, and three times 
that ot Sugarcane. 


(3) A large labor torce can be accomodated and is needed to do both heavy and light 
WOUK. 
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AUTHOK: GU Tingmin [7357 1656 2404) 
GONG Dewen [1362 1795 2429] 
SHEN Jianzhu [3088 1696 2691) 
CHEN Enjiu [7115 1869 0036; deceased] 


ORGS GU and GONG both of the Meteorological Service of Shandong Province; SHEN 
and CHEN both of the Institute of Geography, Chinese Academy of Sciences 


TITLE: "Climatic Regionalization of Drought and Flood in Shandong Province" 


SOURCE: Bei jing DILI XUEBAO [ACTA GEOGRAPHICA SINICA] in Chinese Vol 35 No 3, 
Sep 80 pp 232-24) 


TEX] OF ENGLISH ABSTRAC|: In the past, most of the climatic regionalizations in 
China only showed dry and humid degrees of the atmosphere. In this paper, surtace 
data otf 126 meteorological stations and 90 hydrological stations were used to discuss 
the climatic regionalization of drought and flood in Shandong Province. We intro- 
duce a concept of effective precipitation, utilizing precipitation, potential evap- 
oration and surface runoff to denote a drought and flood index, DW = E/P(i-a). 

Here, E is the potential evaporation, P the precipitation, o the surface runoff 
coetficient and P(l-o) the effective precipitation, According to this index, we 
can get a climatic regionalization of drought and flood but not a climatic region- 
alization of dry air and humidity. 


sinee this index has seasonal variations, we calculated the drought and flood index 
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ot tour di.fterent seasons: 

Spring (March, April, May) DW, = 0,6E/P(l-o) 

Early Summer (june) DW9 = 0.7E/P(1l-o ) 

Mid-Summer (July, August) DW, = E/P(l-o ) 

Autumn (September, October) DW, = 0.7E/P(l-o). 
hese calculating values are used as criteria dividing Shandong Province into 14 
subdivisions. 


ihe index of drought and flood is examined by the water budget and soil and actual 

frought and flood conditions. A table includes DW values of the index of drought 

and tlood compared with actual drought and {lood at Dezou Shi, Lin Xian and Ju Xian 
1963-1972. Classes of drought and flood in the table are taken as the tollowing: 

™ value of 0.5 is a major flood, 0.51-0.75 a flood, 0.76-1.25 normal, 1.26-2.00 a 

frought, and DW value > 2.00 a major drought. We found that the index IW detines the 

limatic regions of drought and flood rather well, 
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AUTHOR: WU Chen [0702 1820) 
MAG 4 Geographical Institute of Hebei Province 


TITTLE: “Relationship of Water and Sand Eruption to Alluvial Sand Bodies in the 
laneshan karthquake Area" 


SOURCE: Beijing DILI XUEBAO [ACTA GEOGRAPHICA SINICA] in Chinese Vol 35 No 3, 
sep 80 pp 251-258 


EX! OF ENGLISH ABSTKACT: During the great 1976 Tangshan earthquake with magnitude 
/.8, a large-scale water and sand eruption occurred on the ground surface of Tang 
han district. Based on field survey and aerial photo interpretation, it is evident 
that the spatial distribution of water and sand eruption, including the eruption 
density, morphological character and various extending orientations and their 
figurations, was associated not only with the sand bodies of ma ro-geomorphic unit 
but also with the micro-geomorphic types, especially the ancient river channels. 
Accordingly, we concluded that the distribution of water and sand eruption is« con 
trolled by alluvial sand bodies and the nature of water and sand eruption is 
letermined by geomorphic types. 


As tor the occurrence of water and sand eruption, in addition to shock force with 
relevant energy, a thicker layer of shallow buried and water-saturated fine silt, 
which was covered or uncovered by clayey loam, is necessary condition. Sine 


[Continuation of DILI KUEBAO Vol 35 No 3, Sep 80 pp 251-2 8) 


one geomorphic 'y pes, such as river flood pialt 2, id and perioheral portions of 
slluvial fans, marine plain and shallow buried ancient river channels, comprise 
the above-described sedimentary structure, and a rong earthquake gives a forcetul 
hock to those geomorphic types, water and sand may easiiy erupt from underground 


(o surface. Therefore, in seismological research, ¢! geomor phologist should 
ccurately demarcate the above-mentioned geomorphic types, whereas the engineers 
should pay particular attention to shock-proof measures on such geomorphic types. 
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AUTHOR: JIANG Hongdi [5592 3163 4574) 
ORL: inatitute of I Wdiation Medicine 


rieues "Unfolding Whole Energy Range Neutron Spectrum from Activation Data by 
Noniineac Programing Method" 


SOURCE: Beijing JITLLANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct ‘oh Dt 24) “JG 


TEXT OF | [SH ABSTRACI: This paper introduces a new method in which nonlinear 
programming is applied to determine the whole energy range neutron spectrum from 
the activation data. According to the different measuring conditions, a set of 
appropriate analytical expressions containing some undetermined parameters is 
selected or developed to deseribe the neutron spectrum to be measured, Then, an 
objective fumetlen ite defined on the basis of the physical requirements of the 
neasur ene | ‘he speetrum and minimized by means of nonlinear programming to ob- 
tain the opt imum selution that satisfies the given reaction rates. Thus, a stable, 
onvetgent and continuous panoramic energy spectrum can be obtained while the 
ippear ance / spectrum oselllation and negative valued fluxes are completely 
ivolLde 
: \ied results of the spectrum measured in channel No 1 of the 
Qiry ven. © ineident neutrons are filtered by natural boron 
t , \NG KUEBAO Vol 1 No 4, Oct 80 pp 241-248) 


7 , } e materia of the foils used are MB» Aly S» Ni, Im, Nay Mn, Co, 











AVIHOR: = None 
ORG Research Group tor JLG Laser Interterometer 


TITLEs "JLG Laser Interferometer for the Automatic Measurement of Highly Accurate 
Gauge Blocks" 


SOURCE: Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA) in Chinese Vol | No 4, 
Oct 80 pp 249-257 


‘FXT OF ENGLISH ABSTRACT: The JLG laser interferometer is a new device for 
automatically measuring the highly accurate gauge blocks. Some advanced techniques 
are adopted in the interferometer, such as stabilized laser, electronic digital 
display, white-light interference fringes for automatically locating the surfaces 
of gauge blocks with a hollow corner cube prism and automatic control, etc. he 
limiting error of measured results for the instrument is 

Biin = t (0,03 + 0.2L) am 
where L is the Length of the measured gauge block in meters. 


AUTHOR: LIU Yunzhu [0691 0336 5053] 
MA Runshan [7456 3387 1472] 
GUO Zhengqiang [6753 2973 1730] 


ORG 3 None 
LITLEs "The Exposure St andard” 


SOURCE: Bei jing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 258-265 


[EXT OF ENGLISH ABSTRACT: An exposure standard apparatus developed in NIM is 
presented. The exposure range of the apparatus is from 0,001 to 10 lux + second, 
the uncertainty of the exposure not exceeding + 0.01 logarithmic unit. It basically 
vtisfies the requirements of sensitometry and can be used to calibrate different 
types of sensit ometers. 





AUTHOK: None 
ORG: Keseacch Group tor Photometric Primary Standard, NIM 
PiThLes "The National Photometric Primary Standard of China" 


SOURCEs Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 266-274 


TEX! OF ENGLISH ABSTRACT: This paper gives a summary of the equipment and experi- 
mental resulte of the photometric primary standard of China. A method has been 
found for reducing the luminance fluctuation of the freezing platinum blackbody. 

The fluctuation may be reduced to within + 0.3 percent from one freezing to another 
one. The luminance of the blackbody at the platinum freezing point agrees with that 
at the melting point within 0.2 percent, Our experience with platinum ingot No 2 
showed that it was not contaminated after being used over 300 times. The total un- 
certainty of the photometric primary standard is t 0.33 percent. 


AUTHOR: CHEN Hongpan [7115 1347 4323] 
HEN Shouren [7115 1343 0088] 
‘HU Zaixiang [5969 6528 4382] 


ORG: All of Harbin Industrial University 


TITLE: "The Evaluation of the Emissivity and Temperature of Cavities at the Gold 


Freezing Point” 


SOURCEs Beijing JLLIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 275-286 


‘EXT OF ENGLISH ABSTRACT: According to the IPTS (68), Planck's formula is used to 
detine the temperature above the freezing point of gold (1064.43°C), At present, 

‘he instrument adopted by most countries to realize the scale is the photoelectric 
pyrometer which is used directly for observing the spectral radiance of the bottom 
of 4 graphite cavity enclosed by pure gold. Because the cepth-radius ratio of the 


cavity is large and the cavity wall is very thin, the error caused by the deviation 
of the apparent emissivity of the bottom and the temperature difference 

tw the internal and external wal!s is negligible. However, in practice, when 
the stable thermal equilibrium is established at the gold freezing point, there is 
alway 1 constant thermal flux radiating outward through the cavity opening which 
causes ‘he inner wall temperature to always be lower than the gold freezing point. 


In the present investigation, the integral equation for the temperature 
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distribution inside the cavity and the relation between the emissivity and local 
temperature have been derived according to the basic ideas of Geist. These equa- 
tions have been solved using a computer for some cylindrical cavities and the 
results are given in this paper. 


AUTHOR: ZHANG Baoyu [1728 2450 0810] 


ORL: None 


ITLEs "The Relations Between the Magnetic Fields of the Helix Current and the 
Current Sheet" 


SOURCE: Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4G, 
Oct B80 pp 287-295 


FX) OF ENGLISH ABSTRACT: By analyzing and comparing the magnetic field distribu- 
tions for both the helix current model and the current model, it has been proved 
that the average of the tield for the helix current over the polar angle is 
exactly equal to the field for the corresponding current sheet, no matter what the 

re of the pitch will be. As a matter of fact, the field for the helix current 

equivalent to the superposition of the field for the current cheet and a small 
inusoidal quantity, which can be neglected when the pitch is small enough. fur- 
vermore, the integral average of the field for the helix current over the 

pherical sample placed at the center of the solenoid is equal to the center field, 

hese results have practical importance for further development of the determina- 
tion of the gyromagnetic ratio of protons in the weak field. 








AUTHOR: ZHANG Zhonghua [1728 6945 5478] 


ORGS None 


PrLITLes "A Matrix Perturbation Method for Calculation of Errors in Precision 
Electromagnetic Measuring Circuits" 


SOURCE: Beijing JTLIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 296-298 


reX? OF CNGLISH ABSTRACT: A matrix perturbation method for the calculation of 
errors in precision electromagnetic measuring circuits caused by small variations 
of cireult parameters is described in this paper. The formulae for the calculation 
are very simple and can be used for finding errors of any order. The convergence 
conditions tor the perturbation procedure are also given and discussed. 


AUTHOR: FENG Xinshan [7458 2450 0810] 
WU Youzhang [0702 1635 3864] 


OR ,? None 
ITLEs "A Broadband Coaxial Low Power Standard” 


fF: Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4 


‘» 
ie t » i pp J9Y9- 402 


EXT OF ENGLISH ABSTKACT: The broadband coaxial calorimeter described in this 

paper is developed as the national coaxial low power standard. fhe calorimeter is 
of the twin-load type, one load being the calorimetric body while the other provides 
1 cetercence of temperature. RF power is measured by comparing with standard DC 


power. It operates in the frequency range of 1 to 8.2 GHz, and the total uncertainty 
+ (0.17™0,26) percent. 


coaxial power meter is calibrated by means of the power equation concept. The 
1libration system i 1 three arm generalized reflectometer, and the total calibrated 
ertointy is better than + 0.27 percent. 











AUTHOR; ZHANG Jueai [4545 0658 2088] 


ORG None 
TITLEs “Some Experiments on the LEC-318 Recommendation for an Artificial Ear" 


SOURCE: Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 303-306 


TEXT OF ENGLISH ABSTRACT; Starting from a review of the IEC-318 recommendation for 
the artificial ear, analytical experimencs were made for seven parameters of the 


analogical network with parameters varying one by one. Experimental curves and 
discussions are given. 


AUTHOR: CHEN Shihua [7115 0013 5478] 
ORG: None 


TITLE: "A Fast Logarithmic Charge Amplifier for Liquid Scintillation Counting 
Device" 


SOURCE: Beijing JILIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 307-315 


TEXT OF ENGLISH ABSTRACT: This paper describes a fast logarithmic charge amplifier. 
Using this amplifier in a liquid scintillation counting device, absolute measure- 
ments of the activity of some radionuclides have been realized. The problem of 
closed-loop stability of logarithmic converters employing a bipolar transistor as 
logarithmic feedback element is discussed. A method for measuring the dynamic 

range of the fast logarithmic charge amplifier is described. Finally, the results 
of absolute measurement of the activity of several radioactive solutions by this 
levice are given. 








AUTHOR: ZHANG Lisheng [1728 0500 3932] 


OKs None 


PiThLes "Vechnical Requirements for VCXO Used in t!\e Cesium Atomic Clock" 


SOURCE: Beijing JLLIANG XUEBAO [ACTA METROLOGICA SINICA] in Chinese Vol 1 No 4, 
Oct 80 pp 316-320 


(EX! OF ENGLISH ABSTRACT: Voltage controlled crystal oscillator (VCXO) plays an 
important role in the atomic frequency standard. The performance of VCXO must be 
excellent so as to maintain the inherent frequency accuracy and stability of the 
atomic standard. Through analysis, technical requirements for VCXO used in the 
cesium atomic clock are put forward in this paper. 
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AUTHOR: HOU Guoben [0186 0948 2609) 
LU Lisheng [4151 7787 3932] 


ORG 3 HOU of Shandong College of Oceanology; LU of the Ministry of Jranspor' 


TITLE: "The Hydraulic Model Study of Vertical Shaft Axial Flow Pumps Used [ot 
Inst allation in Docks" 


SOURCE: Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE O} 
OCEANOLOGY] in Chinese Vol 10 No 3, 15 Sep 80 pp 1-12 


TEX? OF ENGLISH ABSTRACT: Critical values of energy loss parameters are estab- 
lished as a basis for the proportioning approach to channel passages in hydraulic 
models on pump intakes. Models that satisfy the critical limits when adjusted 
experimentally to yield zero swirl of the suction pipe inlet have been found to 
Rive reliable indications of full scale behavior, both in regard to swirl correc- 
tion and freedom from air entrainment. 


[his report suggests general results for use by design engineers engaged in power 
stations, pumping stations or dock design, and illustrates them by several examples 
of successful designs. 


AUTHOR: XI Pan'gen [1153 4149 2704] 
ZHANG Shuzhen [1728 3219 3791] 
FENG Shizuo [7458 1102 460A] 


ORG: All ot the Department of Physical Oceanol ogy 
rITLEs "Mathematical Modeling of Circulation in the East China Sea" 


SOURCE: Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE OF 
UCEANOLOGY ] in Chinese Vol 10 No 3, 15 Sep 80 pp 13-25 


fEXT OF ENGLISH ABSTRACT: Based on the “wind curl-thermohaline gradient equation” 
proposed by Feng, one of the authors of this paper, and taking actual bottom fluc- 
tuation into account, a two-dimensional numerical model of steady-state circulation 
triven by a nonuniform surface wind and the thermohaline effect is presented. The 
model is applied to the East China Sea. The computations of winter circulation 
exhibit considerable agreement with the obvserved current field. The trend of cir- 
culation is almost parallel to isobaths. The Kuroshio exerts great influence on 
the southeast part of the area in question where the swift current flows nearly 
along the isobath for 200 meters, its main part outflowing from the east boundary 
of the region, while the other, the smaller one, flows and branches into the 
[sushima Current and the Yellow Sea Warm Current. In the northern part of the 
waters there exists a counterclockwise “Yellow Sea Gyre System" consisting of the 
northward Yellow Sea Warm Current and the southward China Coastal Current. At the 
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southwest end of the area, the Min-Zhe branch of the Kuroshio together with the 
faiwan Strait Counter Current tlows south in a semi-are pattern along the north 
coast! of Taiwan Island, 


(he influence of boundary forces, the thermohaline effect, wind stress vorticity 
and walter depth on the circulation are estimated respectively. The computations 
stow that the circulation pattern of the East China Sea, under the action of bound- 
ary torces (especially the Kuroshio System), depends basically on the bottom 
topography, with the thermohaline effect playing a determinant role in generation 
of the counterclockwise "Yellow Sea Gyre System." The wind stress, however, only 
icts as an adjustment factor and is of less importance. 


AUTHOR: SHEN Yujiang [3088 5148 3984] 


ORG: Department of Physical Oceanology 
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"Numerical Computation of Tides in the East China Sea" 


SOURCEe: Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE OF 
OCEANOLOGY] in Chinese Vol 10 No 3, 15 Sep 80 pp 26-35 


[EXT OF ENGLISH ABSTRACT: For the East China Sea the hydrodynamical equation 
system ot the co-oscillating tides is numerically solved by giving the values of 
tidal harmonic constants at the boundaries. Based on the computed results the 
o-tidal charts of four main constituents are drawn. With the help of the charts 
‘he propagation behaviors of tidal waves in the region are analyzed, and the 


formation of amphidromic systems and the distribution of tidal characteristic 
quantities ace described. The author is of the opinion that there exists a semi- 
fiurnal degenerate amphidromic system near Jilong, Taiwan. 


omparison between the computations and observations exhibits excellent agree- 
hus suggesting that our results can serve general purposes. 
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‘Kies: Miasile Range Ineatrument Ships in Space 


[SHI [KNOWLEDGE OF SHIPS] in Chinese No 7, 20 Oct 80 


ian intercontinental bailietic missile teat on lt 
vere walting for the splashdown at latitude 7” §, 
luded instrument ships, recovery ships and others, 
jown, the instrument package was ejected, and was 
i diver descending from a helicopter. The role of 


chnology te described on the basis of the ships used 


rogram. These ships, with displacements ranging 
were equipped with S band and C band radare for track- 
wed to describe the procedures in the launch stage, 

} incing stage. 
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ese ship Aramo has attracted great attention because 


ict uts fuel consumption by a quarter. The tw 
may be consumed at the rat« tf 65 tons or more 
efficient screw rotation speed and to use wast« 
in propeller rotating speed require iarger pr 
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m some of the 60 percent tf the total energy 
st. Waste gases at up to 400° C are used ¢t 


r: the steam can be used to drive a steam turbine 
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the Kawasaki KSE system, with energy savings ot 
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AUTHOR: YAO Jiajin [1202 1367 3866] 
ORG: None 
TITLE: “The Tomahawk Ship-Launched Cruise Missile” 


SOURCE: Beijing JIANCHUAN ZHISHI [KNOWLEDGE OF SHIPS) in Chinese No 7, 20 Oct &0 
pp 16, 15 


ABSTRACT: The NATO response to deployment of the Soviet SS-20 is deployment of 
the US SSBS-2 and cruise missiles. The Tomahawk cruise missile can carry a con- 
ventional or 200-megaton nuclear warhead, It has a range of 1,700 miles and high 
accuracy. It can be launched from a submarine, from the Nimitz carrier or fram 
Spruance clase destroyers. Its primary guidance system is Terrain Contou taiche 
ing (TERCOM), in which radar altimeter measurements are compared with a storeca 
digttal map. Near the target, where this system renders the missile more easily 
detectible, an image matching system is used, tripling the accuracy. To defend 
againet thie, the Soviets will have to build many dummy installations, but the 
cruise misatle ts cheap, so many can be launched at targets. In the event of a 
mjor European war, Soviet Backfire bombers will launch AS-6 antiship missiles at 
NALO ships bringing troops and materie. to Europe. The cruise missile is the most 
effective means of destroying the Backfire bases. In addition, because US carriers 
are too few to cover all areas in which the US desires to exert influence, the 
use of crutse missiles aboard destroyers and cruisers in the Indian Ocean will be 
benefictal, 


AUTHOR: LIN Hua [27651 5478) 


RG: None 


riTLi: “The Styx Missile’ 


SOURCE: Beijing JLANCHUAN ZHISH! [KNOWLEDGE OF SHIPS] in Chinese No 7, 20 Oct 80 
pp 22-24 


\BSTRACT: During the 1967 Mideast war, the Israeli destroyer Filat was sun« by 
+ Styx missiles at a distance of 15 miles from Port Said, proving the importance 
this weapon. The Styx, or SS-N-2, is an aerodynamic missile powered by a solid- 
(uel auxtitary motor and a liquid-fuel main rocket motor with a range of 35-40 ka, 
which has been supplied by the Soviet Union to many countries. India used it to 
sink a Paktetanit destrover in the Indo-Pakistan war. It can be mounted on small, 
fast boats such as the Soviet Komar and Osa claes. It flies low (100-300 meters, 
‘tT as low as 10 meters) and is not maneuverable. 8yv the 1973 war, Israel had 
leveloped the Gabriel antiship missile, which sank 9 Arab ships, while none of the 
emili, fast Israeli missile boats was hit. Countermeasures were found to render 
the Styx guidance system ineffective. 
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